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ABSTRACT 


The use of differentially deflected fairwater planes 
to control submarine rolling is studied. Because of counling 
between pitch and roll angles, the snap roll that occurs in 
a high speed hard turn affects the stability of a submarine 
not only in the horizontal plane but also in the vertical 
plane. Direct roll control was achieved by making use of the 
fairwater planes ina differentially deflected mode such that 
they could give counter moment to reduce the snap roll. A 
roll controller was designed as a position and velocity feed- 
back controller. Controlled roll angle improved the depth 


and pitch stability of submarine. 
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RENO TON 


The importance of automatic control of pitch-depth and 
roll of a turning submarine has become more obvious with the 
improvements of high speed nuclear submarines. This study was 
aimed at smoothing the roll of a turning submarine in hiah 
speed, which is known as snap roll and causes the basic problem. 
Because of coupling between pitch and roll, before attempting 
to smooth the snap roll, controlling of pitch and depth in 
high speeds turns was also studied. In this study an auto- 
matic course controller was not considered. 

The methods and techniques, which can be used to reduce 
and smooth the snap roll, can be categorized under two head- 
ings. The first concerns changes in the naval architectural 
characteristics of the designs, such as increasing GM (meta- 
centric height) and reducing sail size. The second categorv 
involves those alternatives which make use of an automatic ship 
control system. Rudder sequencing and speed reduction fall 
under this classification. Reference 1 investigates the affect 
of increasing GM, reducing sail size and speed. In Reference 2, 
Stamps designed on automatic roll controller which makes use of 
rudder sequencing as a function of approach speed and instan- 
taneous roll angles. 

In recent studies, which were made by Naval Ship Research 
and Development Center (NSRDC), it was proposed to enhance the 


control of roll in high speed turns by using the Fairwater 
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Milanes disserentially deflected so that they can be used to 
give a counter moment to reduce the roll. The investigation 

of the differentially deflected Fairwater planes effect was 
part of the project called "Imoroved Control For Advanced 
Submarines". The project was carried out under Program Element 
No. 62754N and Task Area 2F434001. The work unit number was 
1-1563-001-74. In Reference 3, estimating the effectiveness 

of differentiallv deflected sailolanes was investigated as a 
part of project mentioned above. Reference 3 was the inspira- 
tion of this thesis and using differentially deflected sail- 


Planes for direct control of roll was chosen as the desian goal. 
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TH YH EWN TO EFO THE PROBLEM 


Turning characteristics of a surfaced submarine, more 
or less, looks like those of a surface ship. But the situation 
in the submerged position shows big differences. These dif- 
ferences are the result of different naval architectural char- 
actacteristics. Sail structure can be considered the main dif- 
ferences and the main source of roll problems. It is well 
known that when a surface ship goes into a turn, it experiences 
an initial inboard roll. After a very short transient it heels 
outboard. The reversing of the roll angle is primarily due to 
an asymmetric rudder. Since the submarine rudder is generally 
well submerged, the surfaced submarine has an outboard roll 
angle during both the initial and the steadv state phase of a 
turn. For a typical submarine, at moderate speed this outboard 
roll angle is less than 10° and does not cause any big problem. 

If a submerged submarine goes into a high speed turn, 
differences from the surfaced behavior are noticed and the 
problem becomes three-dimensionalized. In the very first phase 
of the turn, the ship has a small lateral velocity v and a 
small rotational velocity r. For a starboard turn, at some 
point along the centerline, Xl, the lateral velocity due to 
r, rXl, equals the lateral velocity - v. (For symbols see 
Figure 1, 2, and Appendix A). The point, at which these two 
velocities are equal to each other, is called the instant 
Cer of turn. At some point, forward of the instant center, 


A», the lateral velocity, -v, gives a velocity component from 
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the starboard side. Along with the ahead velocity component 

of the ship, these velocities contribute an angle of attack 

from starboard. At this point, if the submarine has an append- 
age (sail), a lift force is produced on the sail directed to 
port. Depending on the sail configuration and the large moment 
arm of the sail. This force, directed to the port gives a small 
outboard (port) roll. In many submarines this outboard roll can 
not be felt. 

In the steady state phase the angle of attack shifts to 
the other side (port) and results in an inboard roll. But con- 
trary to the outboard roll this inboard roll is quite signifi- 
cant and at high speeds it can exceed 309 which is considered 
hazardous to both men and eauipments. This inward roll is 
called the SNAP roll. In Figure 5, the roll characteristics 
of the ship to a constant 35° left rudder angle at 24 knots is 
shown. As can be seen, the initial outboard roll (in the 
Simulation, outboard is the starboard side) is very small, 
less than 1°. But the Snap roll reaches SE 

Snap roll is not the only problem that can be faced when 
the submarine goes into a turn. Because of the appendage 
(which is mainly the sail), in the submerged position, sub- 
marines do not have very good hydrodynamics from this body 
structure. When the ship experiences a submerged turn the 
sail mainly destroys the waterflow around the hull. As a 
result a pressure difference is created between the unver 


“Ion of the hull and the bottom portion of the hull, 
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And a normal force results directed downward which bodily 
pushes the ship down. It is recalled that, in a turn, after 
the initial phase, the submarine has an inward roll with the 
deflected rudder. Rolling inward, causes this deflected 
rudder to act like a sternplane which gives the submarine a 
down pitch movement. All of these effects come together and 
cause the ship to rise or dive depending on the peculiaritv 
of its design. As it is stated in Reference 4, "The exact 
mechanism of this is not well understood". In Figures 3, 4, 
and 6, the depth, pitch, and rudder response of the ship, to 


° left rudder angle at 24 knots is shown. It is interest- 


a 35 
ing to indicate here that after a long transient period, the 
submarine starts to rise with the bow down eventhough it dives 
initially. It should be remembered that in this simulation of 
the model, no control surface (stern and fairwater plane) was 
used. The changes of these characteristics will be seen after 
the depth-pitch and roll controller is desiqned. 

mae complexity of controlling the turning submarine 
mainly comes from the coupling between the roll and pitch. 
Since roll and pitch are coupled, to be able to control the 
roll, which is our primary concern, pitch control to some ex- 
tent must be accomplished. In the following sections a pitch 
and depth controller is designed by using only the sternnlanes 


because of the goal of using the fairwater plane as a main 


part Of the direct roll controller. 
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ick COMPUGBR STIMULATION OF SUBMARINE MOTION 


The mathematical model of the totally submerged sub- 
marine is represented by the equations of motion which consist 
of six equations such that each of them represent the force 
and moment equations on and around the three rectangular co- 
ordinate axes. These equations were derived by NSRDC in 
Reference 5 and for convenience was repeated in Appendix A. 

In this thesis, derivation of these equations was not our 
concern. More pronounced knowledge about that can be found in 
References 4, 5, and 6. 

The equations of motion, which represent the motion of 
the totally submerged submarine, in six degree of freedom, con- 
Sist of Hydrodynamic coefficient, in undimensionalized form 


except the m (mass) and -— I I, (three moments of inertia). 


y? 


In Reference 7, Drurey translated these six equations into the 
DSL (digital simulator language) form and originated the com- 
puter program which is the mathematical model of the ship. 

For the necessity of using the hydrodymanic coefficients in 


undimensionalized form, the equations containing I. I T 


E 
both sides are divided by ie and these three moments of 


a! 


inertia were used in undimensionalized form. For the same 
reason the equations containing m were divided by ER 
Since the trim control of the ship was not considered 


in this study the ship was assumed in the trimmed condition. 
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For trimming the ship the hydrodynamic coefficient 
mand M, were set to zero. These two hydrodynamic co- 
efficients represent the force and moment acting on the ship 
when the control surfaces are at zero deflection. 

Since all of the simulations of this study were done 
in a turn, an accurate simulation must include the rudder 
actuator dynamics such that any phase lag could not produce 
any unstability. Actuator dynamics used in this thesis was 
the same as used in Reference 3 with the rudder deflection 
angle 5°/sec., so that at the end of the study it could cive 
acceptable comparison level between the two designs. A block 
diagram of the rudder actuator is shown in Figure 7. 

Since our prime concern is to control the roll as well 
as the depth and the pitch when the submarine is in a turn, a 
series of tests were run at different speeds to validate the 
models originated in Reference 7 by Drurey. The model was 
highly satisfactory as shown in Figures 8 through 25. These 
figures record responses in depth, pitch, roll, yaw, speed 
changes, and rudder angle at 18, 12, and 6 knots. In this 
test, turn was commanded 10 seconds after the simulation 
started, to validate the initial condition response of the 
model. Since the submarine was trimmed by eliminating 2, and 
Me, it is seen that until the turn was commanded all responses 
are zero. Simulation at 24 knots was already shown in Figures 
3 through 6. 

It was observed that, at 6 knots the ship continues 


diving unlike the simulation of 12, 18, and 24 knots. The 
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reason for this was that at low speed enough vressure dif- 
ference to cause the ship to rise slowly could not be created. 
After the validation of the model it was decided to pro- 
ceed with the designing of the controllers. In the following 
sections first, the depth and pitch controller design was stud- 
ied. An optimal control scheme such as in Reference 7 was used. 
The sternplanes was the control surface of the controller, and 


the fairwater plane was saved for use in the roll controller. 
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SS Lee DEPTH AND PITCH CONTROL 


A. OPTIMAL SOLUTION TO THE LINEARIZED SUBMARINE EQUATIONS 

The necessity of controlling the depth and pitch of a 
turning submarine, as was stated before, stems from the coupling 
between the states such that any changes of one of them directly 
affects the other one. For this reason, depth and pitch control 
to some extent must be accomplished so that it can provide a 
good stage for the roll control. There are many ways of con- 
trolling the depth and pitch which depend upon the way of using 
stern and fairwater control surfaces combination. The scheme 
that was used in this thesis was an optimal control way which 
was Originated by Drurevs in Reference 7. To preclude unneces- 
sary repetition, only guidelines of the method and the differ- 
ences of this study are to be discussed here. More pronounced 
knowledge can be found at Reference 7 and 8. 

It can be seen in Appendix A that the equations of motion 
are nonlinear. To be able to use this set of equations, they 
must be linearized. This linearization was done in Reference 7, 
and with one control input (sternplane) they can be written in 


the form of: 


ura /, UZa W y? zas /, 


ze 


se DM  -1-2q 


et 


~mw / 1 I, Mà 
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If they are put in state variable form 


wW nm S P 
- + DS 
A 202 952» Q PD 
Where 
Determ - [ x, (- 20) - meza | 
= NE + ZGMw 
Ar1 | (ty iq) Zw Zéi A Determ 
12 = MUZw + (m-ZW)MW uf, 
12 * Determ 
= = + 7 
A21 [Gy Må) Zq záma | Storm 
= o + = e i " 
A22 E (m-Zw) “a | Dat, YI 
B11 = EXE + zamas | A . Determ 
: : 2 
B12 = [nezas + (m-ZW) nas | u Az TYN 


The problem was thought as a linear regulator problem. The 
general scheme of linear regulator problem is shown in Figure 
26. 


The cost function is 
t 


d 
T= > j Ee: + u?nu ] dt 
Es 
E = X-r and r is the reference vector which is 
ordered depth 0 
EL = = 
ordered pitch 0 


49 





*"uo[qoag 1 
o3vI7nboì 
Xxeourg I 
| Of 19] T[013JU0) 1D0eqpa 
: Əd Ieurad 
I3dO puv que 
ld 


(4) x 


*92 sanbta 





50 





As a result 
E -X can be written. 


Following the Reference 7 and 8 


. -1.T T 

K = -KA+KBR 'B K-O-A K and 
«1 T 

UU = -R B KE can be written. 


If the linear equations are auqmented and but in the state 


equation form with the definition 


Y 


Q 
Xy 0 1 -u 0 X1 0 
X, =10 Ai 0 A51 X + Bii Ds 
X4 0 0 0 1 X, 0 
A 00. A12 IE CH E 
Where 
X4 Depth Rate 
X, = Depth 
X4 Pitch Rate 
X Pitch 


In this study, in order to make the depth rate a state vari- 
able the following transformation was used 
Depth rate = w-u O where w represent the comnonent of 


mem the z direction, = represents the nitch angle. In 


SE 


Reference 7 this transformation was originally 

Depth Rate = w and the differences between these two 
transformations showed itself in the calculation of K values 
which is discussed in the following nages. 


Referring to the cost function 


E 0 0 0 
— -= 0 A O O R= C 
0 0 0 0 
0 0 0 B 
and 
Rai Ron 31 Fax 
m- AA A 
DERE us E 
Kua Boa bar Ras 


The selection of the weiahting factors was a trial and 
error process. The series of weighting factors denending upon 
the relative severity of the influence of the states to each 


other was tested and as a result 


E = Depth Error Weighting = 10 

D = Pitch Error Weighting = 8000 
C = Control Inout Weighting = 100 
A = B = 0 was chosen. 


After solution of n(ntl), differential eguations, 
2 


K values associated with the feedback cain was found. It is 


observed that steadstate values of X are constant which is a 


52 





convenient condition for optimality. Results of this solution 
- are shown in Figures 27 through 36. For comparison the original 
A matrix of Reference 7 with the transformation Depth rate =w 
was also used to calculate K values. One of the results is 
shown in Figure 37. It is observed that it could not reach 
steady state value in the reasonable time period. 
After K values associated with the feedback gains were 

found, necessary gains were found via the optimal low. 

u = -i B K E 

The result was 

“Ue DSAD = -0.31623*Z0ER-1.792*Z0D0T + 

36. 069*PERR = 102.63*P1D0T, 
All of the feedback gains were found at 15 knots speed. 

It is obvious that these gains associated with the hydro- 
dynamic coefficient used are function of the speed. To make 
the gains compatible with the speed range of the ship, they 
must be scaled with the function of the speed. In Reference 2, 
the controller gains adjustment as a function of the speed was 
discussed deeply and found that gains associated with the depth 
rate and pitch role error channel was inverselv proportional to 
the instantaneous speed as gains associated with the depth 
and pitch error channel remains same. After the justification 
of these results 


~ 
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UK 


SET + 102.63*(15), ZPLDOT 
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BoAD = URSS, SDO F 
UK 


36.069*PERR + 1539.45 SPIDOT 


/ UK 
where UK is the instantaneous speed (knots). 
B. AUTOMATIC DEPTH-PITCH CONTROLLER AND SIMULATION OF THE 

RESULTS 

The feedback gains, which are the optimal solution of 
the linearized vertical plane equations must be put into the 
general controller schema in which the submarine dynamics are 
represented in six degrees of freedom (i.e., non-linear 
equations of motion). Because only the simulation of this can 
justify the validity of the solution. The complete modified 
version of the Drurey's depth and pitch controller is shown in 
Fighre 38. As it is shown the sternplane controller (actuator) 
was not designed as a part of the controller and put into the 
controller separately with the plane rate 79 /sec. The dynamics 
of the actuator is the same as Drurey's and Stamps' controller 
used. The same actuator dynamics also was used in the roll 
controller which is to be discussed in the following section. 
The actuator block diagram is shown in Figure 39. To preclude 
stability which occurs as a result of the excessive uses of 
the sternplane and avoid use of the control surface with bia 
deflections the limiter was put into the pitch and depth error 
channels. Pitch error limiter was 10° and depth error limiter 
was 20 feet and found by a trial and error process. 

After the completion of the controller design, various 
tests were run to see the effect of the pitch and depth con- 


troller on depth, pitch, roll and sternplane responses. 


65 









iN IH Vus 


* 13] [013U0) uo3id puy uade 


HOLId 


uITIOULNOD 


Hiddd 


*8€ o4nbr, 





T° USTIONuLNOO| -S+ 
INVIANUALS Ene mita | 0733 oz- 
L ANY 
HLduHGd 7 
*33 07+ 


HOLId 
adá daUo 


Hida 
dudado 


66 





*(J0O3Pn32V SOUeId ASFEMITEY 
O3 SIoJou dq “'*IJ030n32V Seue[d uUIJI93S OL SIoJouj SG) 


*SoUPId A9IJeMITEI] DUY UIO3S 104 Topo xz03e2n3o2y OTTNerpAy *6€ Sanbta 


SAA 





YOUU dao4d 


YO 
T s ayosa 
Ja wo sd 





DAS (43) Suv Ad Ad 
uIad*— a + YO = YO 
of Sd Sd 











67 





The results of tests at 6, 12, 18, and 24 knots to a constant 

35° left rudder angle are shown in Figures 40 through 55. 

When these results are compared with the responses of the ship 

with no pitch and depth controller which is shown in Figures 

3 through 6 at 24 knots and in Figures 8 through 25 at 18, 12, 

and 6 knots. The following comparisons can be drawn: 

1. At 12 and 18 knots controller kept the ship's depth 
and pitch stable after very reasonable transient time. 
The steady state pitch and depth are (5.2)9 - (9.2) feet 
at 12 knots and (ASD - (9) feet at 18 knots. 
2. At 6 knots, the controller failed to keep the pitch and 

depth stable. The ship very slowly kept on losing depth. 
The reason was insufficient speed scalina of the feedback 
gains and needs for another control surface (Fairwater 
planes). Lack of fairwater plane control surface in the 
depth-oitch controller unabled it to keep ship stable at 
low speeds. But in this thesis the primary concern was 
to minimize the snap roll which occurs at high speed. 
In Reference 2 Stamos collected data which shows the neak 
snap roll as a function of the approach speed and ordered 
rudder angle. This data is repeated at the following 
page for convenience (Table I). It is seen that dangerous 
snao roll starts occuring at approach speeds above 12 
knots. For this reason unability of the controller at 
low speed was neglected. Because, a control scheme, which 
at low speed uses Drury's original depth-pitch controller 


(Fairwater plane is part of controller) and at high 
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speed uses this design (Fairwater plane is not a part 

of depth-pitch controller but is a part of the roll con- 
troller) can be designed and switching from one to another 
can be achieved. 

3. At 24 knots the controller again failed to keep the ship 
stable. Depth and pitch response went into oscillation 
with big amplitudes. Original roll resnonse of the ship 
which is shown in Figure 5 was destroyed in the sense of 
decreasing amplitudes of roll oscillation. It did not 
reach stable value and appeared to be oscillating in the 
range of 10°. The reason was the following. The big snap 
roll which was around D initially gave very big dis- 
turbance and oscillation. With these big disturbances 
the controller which uses only the sternolane as a con- 
trol surface was unable to do the job. It was thought 
that, if there had been any controller which could have 
decreased the snap roll (Roll controller) it would have 
been able to stabilize the ship in pitch and depth as 
well as in the roll response. 

Based on the results of these tests, proceeding with the 
roll controller design which would make use of fairwater plane 


as a control surface was decided, 
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V. AUTOMATIC ROLL CONTROLLER 


In the introduction of the thesis, the possible choices 

of controlling roll was-discussed briefly and stated that 
fairwater planes were to be used as control fins to control 
the roll. The reason that brought uv this idea lies behind 
the needs of more control surface than was available in the 
present design criteria. Roll control of a high speed sub- 
marine has continued to be a problem over the years because 
of limited control surface. Control engineers have been re- 
stricted using rudder, stern and fairwater plane. All of the 
control surface was not meant to control roll. In the present 
desian criterias stern and fairwater planes are meant to con- 
Brol NE rine in the vertical plane motion (depth-oitch). 
The rudder has been used as a part of course controlling of 
the ships. Among the three control surfaces, rudder is the 
only one that has direct correlation with the submarines roll 
angles. In the latest study, which was done by Stamos 
/Reference 2/7, it was used to control roll. The concept of 
Stamps roll controller was based on the idea of adjusting the 
initial rudder angle order as a function of the integral of 
error between the allowed maximum roll and actual roll anales. 
The initial rudder order was chosen such that the peak roll 
expected for a given aporoach speed would be less than the 
maximum allowed roll. The integral of roll error was then 


computed and scaled to represent and additive term applied to 
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the initial rudder order. Stamps design didn't give any 
structural changes to the present navy submarines design 
criteria. In a sense of simplicity it was perfect. But it 
prohibited using hard and excessive rudder order which is 
highly desirable at some submarines reauired maneuvering in 
certain tactical areas. Since the Stamos design slowed down 
the yaw rate and caused the ship to change its course very 
slowly, alternative design ways were investigated and using 
fairwater plane in the differentially deflected mode was 
chosen as a possible improvement. 

The concept of roll control by means of fairwater planes 
is based on the idea of deflecting the fairwater vlane dif- 
ferentially such that it can give a roll moment in opposite direc- 
tion to the instantaneous roll angle. If the ship has a 
roll angle to starboard side, the planes are to be deflected 
to give a roll moment to portside. The positive sense of this 
additional roll moment created by the differentially deflected 
Sailplanes is the same as of Peference 5. The positive sense 
of sailplane deflection angle adopted in this thesis (positive 
for port sailplane deflected leading edge up, starboard sail- 
plane deflected leading edge-down) is in agreement with that 
used for other control surfaces in Reference 5: Positive 
when the surfaces are deflected in such a direction as to in- 
Crease the relevant angle of the submarine about its mass 
center. In the case of differentially deflected sailplanes, 
the angle is the roll angle and this is defined in Reference 5 


to be positive starboard side down. 
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Since, in the present design criteria of navy submarines, 
sailplanes are not used in a differentially deflected men 
standard eguations of motion which are developed by NSRDC in 
Reference 5, do not have correlative terms which take the 
additional moment term, created by the sailplanes into account. 
For this reason, before starting to design the roll controller, 
this additive moment must be estimated. Once the counter 
moment is found it can be placed in the righthand side of the 
roll moment eguation. 

A. ESTIMATING ROLL MOMENT DUE TO DIFFERENTIALLY DEFLECTED 

SAILPLANES 

In vertical plane motion, the normal force due to the 


fairwater planes (sailplanes) is given by the equation of 


ML e E 


Where 9 = 
L = ship length in feet 
U = forward speed 


26h : hydrodynamic coefficient associated with the 
Sailplane deflection (in conventional mode - right and left 
Side moves together in the same direction). 

6b = Deflection of the sailplane in radians. 

After setting =2 the normal force equation becomes 

Fy = 12.02.25, 0b 


This force is due to both sides of the sailplanes. The 


force due to one side of the plane pair is half of the total 





force that is 5 L -ut zog OR If both sides of the planes 


are deflected by the same amount and in opposite direction 
compared to each other, this configuration creates a moment 
associated with the moment arm between these opposite forces. 
This is shown in Figure 56. Total amount of the moment due to 


sailplanes can be written 


M H Ig 5b. (moment arm) 


A (moment arm) *É 


A 
il 


Lin, Ob» (255 *moment arm), 


if this momemt term is placed in the right side of the 
equation of motion about the body axis system X-Axis (roll 


axis), the complete equation takes the form of 
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Mm this last form Opis replaced by b. to indicate that the 
fairwater plane is to be used in differentially deflected mode 
( Sprepresents deflection anale of fairwater plane in con- 
ventional usage. Both right and left sides move in the same 


direction), 
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bagure 56.1, Positive Direction Of Differentially 
Deflected Sailplanes. 


.2. Negative Direction Of Rolling Moment. 
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The moment arm used in this study was estimated by in- 
voking the data from Reference 6 that belongs to the un- 
Classified fictitious submarine. In the reasonable proportion 
to the length of the ship, the moment arm was chosen as 7 ft. 
and, according the data, it was thought acceptable. 

If the new hydrodynamic coefficient is defined associated 


with the fairwater planes in differentially deflected mode 


2C p _ (ZSto'*moment arm) 


L 
“SF = (0.00558 x re 
= 0.0035517. 


After the estimation of ZC pe as a hydrodynamic co- 
efficient, it is placed in the modified equation discussed 
above. After these modifications, in the equations of motion 
Ot the term which is a function of the Öb (sailplane des 
flection in the conventional mode) was set to zero by defining 
Sb= 0. As result, the equation of motion about the X-axis 
(roll axis) has correlative term between the fairvater plane 
and roll angle of the ship so that this additive moment term 


can be used to give feedback to smooth the snap roll. 


B. ROLL CONTROLLER DESIGN 

In contrast to the depth and pitch controller, the roll 
controller was designed by using the nonlinear equations of 
motion and linearizing was not attempted. In the orevious 
section, the modification of the roll moment equation due to 


the imoosition of moment term which stems from the fairwater 
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planes deflection was discussed. All through the roll con- 
troller design this modified version of the equations was 
used. The reason that led to using six nonlinear equations 
was such that linearizing of this equations in six dereee of 
freedom was found difficult because of the terms coupled with 
each other and mainly was thought that from the output of the 
system dynamics the roll angle and the roll rate (p) could 
give enough feedback information to accomolish the compensation. 
The principles for designing the roll controller were 
Simple. If the shiv has a positive roll angle (starboard) 
the fairwater plane is to be deflected so that it can give 
negative roll moment, i.e., starboard sailplane is to be de- 
flected leading edge up while oort sailplane is deflected 
leading edge down or vice versa. As long as the ship has a 
roll angle this would cause the fairwater plane deflection by 
way of the feedback channel used. In the design it was assumed 
that, the fairwater planes actuater was capable of givina e- 
gual amounts of deflection command in opposite directions to 
both starboard and port sailplanes. Since the design was to 
invoke extensive simulation study (since nonlinear equations 
were used) in determination of the feedback parameters, the 
first basic controller attemot was the proportional controller 
because of its simplicity to design and implement. The pro- 
portional controller is shown in Figure 57.1. Referring to 
the figure, if the submarine has a roll angle this would give 


the position feedback to force the system to reach zero roll 
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Figure 57.1. Position Feedback Controller. 


.2. Position And Velocity Feedback Controller. 





angle position. Determination of Kl (proportional constant) 
was a trial and error process. Even though position feed- 
back alone gave big improvement in the system dynamics, as is 
to be discussed in the following pages, it failed to stabil- 
ize the system in some overating conditions that are consider- 
ed very likely to be faced. This leads to the compensation 
of the system by velocity feedback and other modifications 
(Limiter in the position feedback channel). In the following 
pages the design procedure which uses extensive computer sim- 
ulation is discussed from the simplest case of proportional 
controller to the last modification that TL Ed the system 
in various operating conditions and over a wide range of speeds. 
l. Proportional Gente lier 

The first attempt to control the roll was to design 
a proportional controller as is shown in Figure 57.1. Refer- 
ring to the figure, it is seen that ordered fairwater plane 
deflection m function of the roll angle and the proportional 
constant Kl. Such that 

ORDERED FAIRWATER DEFLECTION = DFOD = -K1*ROLL 

Since the reference signal is zero the system always looks for 
zero roll angle. Ordered fairwater deflection as an input to 
the fairwater planes actuator causes deflected planes at the 
output due to the actuator dynamics. Actuator dynamics were 
the same as those of the sternplanes actuator shown in Figure 
39. With the counter roll moment created differentially de- 
flected fairwater planes, system dynamics try to decrease the 


roll angle. 
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Determination of Kl was a trial and error procsss. 
The stability range was 1 through 4. By inspection of the 
results of computer simulations, Kl = 3 was chosen as a best 
choice in the sense of snap roll and steady state value of 
the roll responses. With the determined Kl = 3 value the 
system was tested at 24 knots approach sneed to a 35° rudder 
command and the results are shown in Figures 58 through 62. 
Before analyzing the results it is instructive to indicate 
here that all through the design procedure the worst condition 
that could happen was always taken into account. For this 
reason the controller was designed at 24 knots base sneed and 
to a 35° constant rudder angle which could give the worst snap 
roll in the range of sneed of interest. Bv lookinc at the 


Figures 58 through 62, the following results can be summarized: 
1. As was predicted, the snap roll decreased to almost 


9.6 from 37? and roll response reached a steady state 


4 
value of 39.5, In Figures 52, 53, and 54 it was shown 
that before implementation of the roll controller the 
system was unstable to an identical test. Adding the roll 
controller made the system stable at hich speed. In depth 
and pitch responses great imorovement has been made. 
Steady state values of the depth and pitch were almost 

8 feet and 49.5 which were acceptable. While this deoth, 
pitch and roll control improvement has been obtained, the 


stern and sailplane was used moderately such that they 


didn't reach any saturation (+ d 
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Figure 61. Sternplane Angle vs. Time. Kl = 3. UCK = 
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2. Investigation of the roll response in Figure 60 lead 
to further study. In the steady state value of the roll, 
the oscillation with almost fixed amplitude and fre- 
quency was observed. The amplitude was around 041. As 
a first guess, the oscillation was thought to be a stable 
limitcycle. To expose the problem the system was tested 
with initial roll response a) inboard roll angle) at 24 
knots to a 35° constant rudder angle. Figures 63, 64, 
and 65 record the depth, pitch, and roll responses. The 
system became unstable. The svstem failed to compensate 
itself when it was commanded to turn with the 35° rudder 
deflection and 5° initial inboard roll. Any submarine 
maneuvering in the tactical area may experience a turn 
with initial roll angle. Disturbances from the heavy sea 
state might cause the submarine to roll eventhough it is 
in a straight course. If the submarine is commanded to 
turn at this moment, it was shown that its control is lost. 
That is why the situation was thought unacceptable and com- 
pensation of the system was attempted. It was thought that 
the reason for the failure was lack of enough feedback in- 
formation and compensation of the system with velocity 
feedback was attempted. 
The following section describes the modification and improve- 
ment obtained. 
2. Compensation Of The System With Velocity Feedback 
The modified controller block diagram is shown in 


cure 57.2., and DFOD = -Kl1*ROLL-K2*P, 
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It was felt that the velocity feedback should give to the 
system better dampina and better stability characteristics. 
Determination of the K2 value was again a trial and error 
process and the stability range was found to be 1 to 15 by 
computer simulation. Bv inspection of the results of the 

tests with different K2 values, K2 = 10 was considered the 

best choice in the sense of snap roll (max roll) and steady 
state value of the roll responses. Since the reason for intro- 
dncina velocity feedback was instability in the presence 

otf an initial roll angle, tests with exaggerated initial roll 
(209 inboard) were made. It was previously denoted that the 
ship has almost Dos snap roll at 24 knots to a 35? rudder 
command. To force the ship to turn with a hard rudder command 
when it already has a very serious initial roll angle was 
thought a good example of the capability to control the ship 
in three dimensions. In fact, any submarine cruising in a 
Stratgnt course is unlikely to face this much roll angle from 
heavy seastate. In Figure 66 the result of this test is shown. 
Before the implementation of the velocity feedback the system 
was unstable to a turn command with initial 5° roll angle. 


0 initial inboard roll. The 


And now it is stable even with 20 
result was such that steady state roll angle was almost 49.5 
and before reaching steady value on oscillation with decreasing 
amplitude was observed. In all of the aforementioned tests, 
the ship response was investigated to a steady 35? rudder 


command, The rudder was commanded to a ES deflection at the 


beginning of the simulation and after reachina 359 it was not 
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changed to any other command. But any submarine maneuvering 
in tactical areas can experience successive rudder commands 

in opposite directions. For this reason the test, in which 
the rudder deflection was commanded to 35° in the time inter- 
val of 20 to 80 seconds was required. Figures 67 through 71 
record the depth, pitch, roll, sternnlane and sailvlane de- 
flection. In this test the initial roll angle was again 20° 
inboard. The result was unstable. Inspection of these 
figures reveals a very important reason of the failure. Curve 
number 2 in these figures represents the rudder response. But 
Since in es vertical axis the automatic scaling was used 
associated with the output responses of interest, in some 
figures the rudder response in the time interval of 20 to 80 
seconds was seen less than 359. For this reason, to overcome 
misinterpretation, curve number 2 should be interpreted as the 
time interval where rudder deflection was 35°, Bv inspection 
of the figures 67 and 68 its seen that the depth and pitch re- 
sponse are almost unchaged. The very small MEG are due to 
the roll oscilation which stems from the initial roll angle. 
In Figure 70 it is shown that in the time interval of 0 - 20 
seconds the sternplanes oscillate with very small amolitudes 
to compensate these depth and pitch changes. But in the same 
time interval (0 - 20 seconds) the sailplanes oscillate with 
big amplitudes to overcome the roll oscillation started with 
the initial 20° inboard roll angle (Figure 71). Because of 


the roll controller, the initial roll angle has more effect 
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on the sailplanes than on the sternplanes. In the time inter- 
val of 0 - 20 seconds of Fiaure 69 the roll oscillation started 
with the initial roll angle is seen. Since the direct effect 
of the roll angle on the sailplane deflection, this initial 
roll oscillation causes the sailplane to oscillate and this 
oscillation causes initial roll oscillation not to be damned 
out. As this is going on, at time equal to 20 seconds the 
rudder was commanded to 35? full deflection and this gives more 
unstability. After analyzing this and previous results the 
following conclusion has been reached. Since the roll angle 
was directly fedback via Kl, the sudden and big roll resnonse 
changes cause the sailplane to oscillate with big amplitudes 
and this leads to unstability. To overcome this difficulty a 
limiter was placed in the roll error feedback channel with the 
magnitude of +5 (The same reason had lead to the placing 
depth and pitch error limiter in the depth and pitch controller 
design as was discussed in Section IV.B.). After this last 
modification the system was tested under various conditions 
which are to be discussed below. The complete controller block 
diagram with the depth and pitch controller is shown in Figure 
72. By inspection, the simulation results of the last version 
of the controller can be sumarized as follows: 
l. In Figures 73 through 87 the depth, pitch, roll, stern- 
plane and sailplanes responses at 24, 18, and 12 knots to 
35° rudder command in the time interval of 20 - 80 seconds) 


with the initial 20° inboard roll angle are shown. The 
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same test without roll limiter as was discussed akove 
gave an unstable system but after inserting the limiter 
the result was stable over the speed range of interest. 

In contrast to the results shown in Figures 67 through 71, 
there is no oscillation in the roll and sailplanes re- 
sponses in the time interval of 0 - 20 seconds and it 

was observed that the limiter gave very good damping to 
the system. This last statement can be justified by 
inspection of Figure 75 such that at time equal 20 

(when the rudder is commanded to full deflection) the 

roll response is almost zero without any previous oscil- 
lation. As the stabilization in three dimensions was being 
reached the stern and sailplanes never went into saturation. 
Figures 88 through 102 record the depth, pitch, roll, 
sternplane and sailplanes responses at 24, 18, and 12 
knots to 35° rudder command (the rudder is commanded to 
full deflection of the beginning of the simulation) with 
zero initial roll angle. It is ENE that without 
using excessive stern and sailplane deflection areat sta- 
bility in the three dimension has been reached. As was 
shown before, Figures 44 through 55 record the responses 
of the system without roll controller to the identical 
test. Comparison of these tests is made in the Table II. 
It is seen that at high speed great improvement in the 


sense of maximum and steady state value has been obtained. 
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3. Figures 103 through 105 record the roll responses of the 
system at Se, and 12 knots to 35° rudder deflection 
(the turn was commanded at the beginning of the simu- 
lation) with 20° inboard roll angle. The results are 
stable. Comparison of Figure 103 with Figure 66 shows 
the improvement obtained by inserting a roll limiter such 
that it damped out the oscillation. 

Tt was shown that utilizing the fairwater planes as 
the control surface of the roll controller not only stabilized 
the roll response but also qave big improvement in the depth 
and pitch controller. In the following section stability 


tests of the svstem to disturbance moments is discussed. 


Se ROME cONCROLEER STABILITY TESTS 

Since the roll controller design was carried out usina 
nonlinear equations in six degrees of freedom, the stability 
analvsis of the system was not practical. To investigate the 
stability computer simulation was used by applying disturbance 
moments. The magnitude of the moments was 

DISTURBANCE = KS* (UCK) “*L? where 
KS = 4.0*0.0003*576/ (UCK)? 

The disturbance moment was calculated by assigning non- 
zero value of KS which was normally zero and represents the 
hydrodynamic coefficient which gives the rolling moment when 
body angle and control surface angles are zero. The magnitude 


of the disturbance was four times bigger than that used in the 


stability analysis of Stamps controller. At 24 knots the 
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magnitude of the disturbance was 1. 1028x10” Pete such that 


the Stamos controller was unable to stabilize the system to 

the disturbances of the magnituce in the order of 10’. Figures 
106 through 120 record the depth, pitch, roll, sternplane, and 
Sailplanes responses to the inboard roll moment dis- 

turbance mentioned above. The distrubance was applied at time 
equal 20 as a step moment. Curve number 2 represents the dis- 
turbances. Inspection of the results shows that the system is 
Stable in three dimensions to the disturbances moment. The 
worst affect of the disturbances was seen at low speed test 
which was normal. But even in this case, after a transient re- 


soonse the system was stable. 
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Figure 106.1. Depth vs. Time. Roll Stability Step Test. 
UCK = 24 Knots. Rudder Ordered = a 


.2. Disturbance Moment vs. Time (1.1028*10/ 


ft-lb. ) e 
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Figure 108.1. Roll vs. Time. Roll Stability Step Test. UCK 


= 24 Knots. Rudder Ordered = 35. 
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.2. Disturbance Moment vs. Time (1.1028x10 ft-lb.). 


156 


C 


eu. 


X SCALE = 40 Seconds Per Inch. 
Y SCALE (1) = 4 Degrees Per Inch. 


00 4.00 8.00 FO 16.00 20.00 


0» 


4.00 


nuc DESEE 30.00 Ia 160.00 200.00 


Figure 109.1. Sternplane Angle vs. Time. Roll Stability Sten 
Test. UCK = 24 Knots. Rudder Ordered = 35°. 


.2. Disturbance Moment vs. Time pa re 
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Figure 110.1. Sailplane Angle vs. Time. Roll Stability Step 
SS ——'— 24 Knots. Rudder Ordered = 35% 
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Esocure 111.1. Depth vs. Time. Roll Stability Step Test. 


UCK = 18 Knots. Rudder Ordered = 35°, 
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ure 112.1. Pitch vs. Time. Roll Stability Step Test. 
UCK = 18 Knots. Rudder Ordered = 35°, 
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ure 113.1. Roll vs. Time. Roll Stability Step Test. 
UCK = 18 Knots. Rudder Ordered = 35°, 
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Figure 114.1. Sternplane Anale vs. Time. Roll Stability Sten 
Test. UCK = 18 Knots. Rudder Ordered = 35°. 


.2. Disturbance Moment vs. Time. (1.1028x10/ft-lb.). 
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Fioure 115.1. Sailplane Angle vs. Time. Roll Stability Steo 
Test. UCK = 24 Knots. Rudder Ordered = Sov. 


EE vs. Time. (1.1028x10/ ft-lb.). 
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Figure 116.1. Depth vs. Time. Roll Stability Step Test. 
UCK = 12 Knots. Rudder Ordered = 35°. 
.2. Disturbance Moment vs. Time. (COE ACH 
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Mies 117.1. Pitch vs. Time. Roll Stability Step Test. 
UCK = 12 Knots. Rudder Ordered = 35°. 


2. Disturbance Moment vs. Time. (1.1028x10/ 


ft-lb.). 
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Figure 119.1. Sternplane Angle vs. Time. Roll Stabilitv 
Step Test. UCK = 6 Knots. Rudder Ordered = 35°. 


‚2. Disturbance Moment vs. Time. (1.1028x10'*t-1b.). 
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Figure 120.1. Sailplane Angle vs. Time. Roll Stability Sten 


Test. UCK = 12 Knots. Rudder Ordered = Qo 
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VI. DISCUSSION, CONCLUSIONS AND RECOMMENDATIONS 
LO NEW DHER WORK 


A. DISCUSSION AND CONCLUSION 

The characteristic of a submarine which is totally 
submerged and in a high speed turn was discussed in Chapter II 
and the simulation results were given in Chapter III. The 
unhydrodynamic hody structure when it is submerged, which is 
mainly due to the appendades (sail), was considered the main 
source of the problem. The necessity of the sail that provides 
rooms for periscove and other vital controls that are located 
in the tower does not vermit a small sail structure. On this 
acount, a control scheme must be used to compensate the prob- 
lems that stem from the sail. The complexity of controlling a 
submarine in a high speed turn, comes from the coupling between 
states and three-dimension dynamics. In Chapter IV an automatic 
depth and pitch controller was presented briefly with the modi- 
fication of the original design represented in Reference 7. 
The automatic control scheme gave the advantage of taking depth 
WD Ech error and their rate into account. But, since only 
the sternplane was used as a control surface, the results made 
the roll controller design necessary for controlling the shin 
within the wide range of speeds and under various disturbance 
conditions. The imposition of the fairwater plane to control 
roll, which was the main goal of this thesis, qave a great 


amount of expected improvement to the system. The roll 
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controller scheme was carried out by using equations of motion 
(nonlinear) in which the vertical plane dynamics were stabil- 
ized by the imposition of the depth-pitch controller, and roll 
state and roll rate gave enough information to provide closed- 
loop control for the ship. The control scheme was a direct 
roll controller compared to the controller desianed by Stamps 
in Reference 2. Stamps's controller was an indirect roll con- 
troller and in his design the rudder was used as a function of 
the integral of roll error which was the difference between the 
allowed maximum roll and the actual instantaneous roll angles. 
The initial rudder ordered was to be such that the peak roll 
expected for a given approach speed would be less than the 
maximum allowed roll. The integral of roll error was then to 
be comouted and scaled to represent an additive tern applied 
to the initial rudder order (8ro). The instantaneous rudder 


order (DRCOM) was then to be computed 


DRCOM = ro + Ky O o) at: 


The technique reported in this thesis has the benefits of 
simplicity and costs much less than alterations of the desion 
criteria of Navy submarines. But since it prohibited use of 

a hard rudder deflection (35°) in the early phase of turn, the 
time required to pass any specified yaw angle was quite long. 
For this reason, the roll control, ina direct sense, by making 
use of fairwater planes as control fins, was investigated in 
the proposed design. As a result, not only were the dangerous 


snap rolls decreased to magnitudes of ES but also with the hard 
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rudder deflection a big van rate was achieved. In Table III, 
the verformance of the submarine with each of these two roll 
controllers is compared. In the Stamps' design, maximum roll 
angle is oredetermined and the controller is designed not to 
exceed this value. The higher allowed maximum roll angle gives 
higher vaw rates. In the proposed design, maximum roll angle 


is the result of the system's dvnamics. 


TABLE III 


COMPARISON OF THE PROPOSED AND STAMPS'DESIGN 


NEE, Ob MAX APPROACH TIME REQUIRED IN SECONDS TO 
CONTROLLER ROLL SPEED PASS SPECIFIED YAW ANGLE 








STAMPS' DESIGN 24 52 91 7100 
12 39 67 7100 
12 27 47 83 
O 
BE un 55.61 24 Nas 20 41 
FATRWATER 4~.69 18 14 25 54 
PLANE Bo A 12 19 BG 80 


It is seen that the use of differential fairwater planes 
provides a faster response with reduced maximum roll anale at 
all sveeds. 

In the proposed design, only one fairwater nlane actu- 
ator was used and it was assumed that it gives deflection in a 
differential mode. Switching criteria from the conventional 
usage of the fairwater plane (both starboard and port sides 
move simultaneously in the same direction) to the differential 


mode was not discussed, 
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This research has shown that imposition of the fairwater 


Enoc as a part of the roll controller improved the control of 


a turning high speed submarine in three dimensions. Dangerous 


O : 
snap roll was decreased to around 5 and with only sternplane 


and the stabilized roll response high improvement in the depth 


and pitch characteristic was achieved. 


B 


RECOMMENDATIONS FOR FURTHER WORK 


The following research would be worthwhile in future 


Studies: 


l, 


Automatic depth and pitch controller by using only 
sternplane in the sense of an optimal feedback controller 
for linear tracking should be designed and the combina- 
tion of this design with the present roll controller 
should be investigated. This study would give good in- 
Sight into the problem which comes from the depth changing 
when the submarine is in a high speed turn. 

Switchina criteria from the roll controller of Reference 

2 to the proposed roll controller should be established 
and these two designs should combine together. This would 
give better depth and pitch control to the ship when it 
has low speed where a great amount of roll control is 

not needed. 

With the coordination of NSRDC a more accurate hydrodynamic 
coefficient associated with the differentially deflected 


Sailplane should be obtained. 
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ma eron Of an integral of error to the roll control 
should be studied to see if the steady state roll angle 


can be reduced. 


ig 


_ APPENDIX A 
EQUATIONS OF MOTION 


The following set of equattons are referred to a body 
fixed system of axes which are cotnetdene with the prin- 
cipal axes of fnertía of the body. The origin of this 
axís-syStcm is located at the assumed center of mass of 
the bodv 


Equation of Motion Along che Body Axis System x-Axís 


m(u - vr + wq)= 5 (7 LX aq air E ' rp] 
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Equation of Motion Along the Body Axis System y-Axts 

m(- wp*ur e Zo Te: Pe Yat] 
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mid - uq * vp) - 7 


Equation of Motion 
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Along the Body Axis System z-Axíis 
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Equation of Motlon About the Body Axis System x-Axis 
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Equation of Motion About the Body Axis System y-Axis 
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Equation of Motfon About the Body Axis System z-Axís 
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AUXILIARY EQUATIONS 
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NOMENCLATURE 


All symbols used in the equations of motion and ín the 


auxitlfary ecuations and relationships which appear Ín this 


report are defined below. 


Any dimensions involved will be 


consistent with the foot-ponnd-second system oC units. All 


angles are in degrees. The Fortran varíables correspondíng 


to these symbols are shown in Appendix B . 


SYMBOL 
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DEFINITION 


A dot over any symbol signi- 
fies differentiation with res- 
pect to time. 


Buoyancy force which is posi- 


tive upwards. 


Mass of the submarine tncluding 
the water in the free fJ5odfag 


spaces. 
Overall length of the submarine 


Lincar velocity of origin of 
bady axes relative to an earth- 


fixed axis systen. 


Component of U along the body 


x-axis. 


Component of U along the body 


y-axis. 
e 


Component of U along the body 


z-axis. 


Ve 
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Command speed: A steady value 
of u for a given propeller rpm 
when A, B and control surface 

angles are zero. Sign changes 


with propeller reversal. 


Longitudinal axis of the body 


fixed coordinate axis system. 


Transverse axis of the body 


fixed coordinate axis syste. 


Vertical axis of the body fixed 


coordinate axis systen. 


Distance along the x, axis of an 


earth-fíxed axís systen. 


Distance along che y axis of an 


earth-fíxed axís system. 


Distance along the d axís of an 


earth-fixed axis system. 


Component of angular veloctty 


about the body fixed x-axis. 


Component of angular velocity 


about the body fixed y-axis. 


Component of angular velocity 


about the body fixed z-axís. 


The z coordinate of the center 
of buoyance (CB) of the subma- 


rine. 





9b, Ds 


ôr , Dr 


35 , Ds 


CR -O. ea o 


x 
re 


tí 


Angle of attack, 
Angle of deife, 


Deflectton of bowplane ( or 
eailplane ) 


Deflection of rudder. 
Deflection of sternplane. 

The ratio u./u. 

Pitch angle. 

Yaw angle. 

Roll angle. 

Mass density of sea water, 
Weight of water blown from a 
particular ballast tank Ident- 
ified by the integer assigned 
to the Index t. 

Angular velocity. 


Time. 


Location along the body x-axts 


of thc center of mass of the (th 


balÎast tank uhen this tank is 


filled with sea Water, 


Er 


. SA)" a 


L 
Rye 


(ES) 


Sp 


/ i 
Ky M , A uw Ks. 
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Propulsion force (see auxt- 


liary equations and relation- 


ships). 


Monent 


marine 


Moment 


maríne 


Monent 


narine 


of inertia of a sub- 


about the x-axís. 


of tnertía of a sub- 


about the y-axis. 


of inertia of a gub- 


about the z-axis. 


Non-dimenslonal constants cach 


of which ts assigned to a parti- 


culer force term in the equation 


of motion about the body x-axis. 


Non-dimensfonal constants each 


of which ts assigned to a parti- 


cular force term tn the equation 


of motion about the body y-axis. 


Non-cinenstonal constants each 
of which ís asaigned to a parti- 


cular force term In the equation 


0 * $ » 
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uu ' “vv ° “ww ° rde Y “Sate ” 
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d ' d ' cular force term in the equation 
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